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Field-emission patterns from a clean tungsten tip apex induced by femtosecond laser pulses have 
been investigated. Strongly asymmetric field-emission intensity distributions are observed depending 
on three parameters: (1) the polarization of the light, (2) the azimuthal and (3) the polar orientation 
of the tip apex relative to the laser incidence direction. In effect, we have realized an ultrafast 
pulsed field-emission source with site selectivity of a few tens of nanometers. Simulations of local 
fields on the tip apex and of electron emission patterns based on photo-excited nonequilibrium 
electron distributions explain our observations quantitatively. Electron emission processes are found 
to depend on laser power and tip voltage. At relatively low laser power and high tip voltage, 
field-emission after two-photon photo-excitation is the dominant process. At relatively low laser 
power and low tip voltage, photoemission processes are dominant. As the laser power increases, 
photoemission from the tip shank becomes noticeable. 

PACS numbers: 79.70.+q, 73.20.Mf, 78.47.J-, 78.67.Bf 



I. INTRODUCTION 

Field emission from metallic tips with nanometer 
sharpness has been introduced some time ago as highly 
bright and coherent electron source Only recently, 

pulsed electron sources with high spatio-temporal reso- 
lution were realized by laser-induced field emission from 
such tips (9l-[T2|. Potentially, spatio-temporal resolution 
down to the single atom and the attosecond range ap- 
pears to be possible 0, [H, [lo|- Such electron sources 
would be very attractive for applications in time-resolved 
electron microscopy or scanning probe microscopy. How- 
ever, the interaction of the laser pulses with the sharp 
tip and the electron emission mechanism are not yet fully 
understood 

When a focused laser pulse illuminates the tip, optical 
electric fields are modified at the tip apex due to the exci- 
tation of surface electromagnetic (EM) waves that couple 
with collective surface charge excitations to form e.g. sur- 
face plasmon polar itons. Interference effects of the result- 
ing surface EM waves can lead to local field enhancement 
[l^ . Depending on the field strength, different electron 
emission processes become dominant [9]. For relatively 
weak fields, single electron excitations by single- or multi- 
photon absorption are dominant, and photo-excited elec- 
trons are tunneling through the surface potential barrier; 
such processes are termed photo-field emission. On the 
other hand, very strong local fields largely modify the 
tunneling barrier and prompt the field emission directly, 
leading to optical field emission. So far, the different 
emission processes were disputed in the literature, while 
the local field enhancement was treated as a static effect 
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such as the lightening rod effect [ol, [To|, fT5l - [T7| . Hence, 
local fields on the tip apex are considered to be symmet- 
ric with respect to the tip axis. However, when the tip 
size is larger than approximately a quarter wavelength, 
dynamical effects are predicted to occur [18]. 

Here, we used a tip whose apex was approximately 
a quarter wavelength and we investigated laser-induced 
field-emission patterns. We have found that dynami- 
cal effects substantially influence the symmetries of lo- 
cal field distributions and thereby field-emission intensity 
distributions. Varying the following three parameters 
changes these distributions substantially: (1) the laser 
polarization, (2) the azimuthal and (3) the polar orienta- 
tion of tip apex relative to the laser incidence direction. 
These are effects that had not been observed in earlier 
experiments [i, Hoj. At the same time, we realized an 
ultrafast pulsed field-emission source with emission site 
selectivity on the scale of a few tens of nanometers. In 
our previous paper, simulations confirm that the photo- 
field emission process is dominant in laser-induced field 
emission [20|. Here, we further discuss electron emis- 
sion processes and their dependence on laser power and 
tip voltage by investigating electron emission patterns, 
Fowler-Nordheim plots, and calculated electron energy 
distributions. At relatively low laser power and high tip 
voltage, field emission after two-photon photo-excitation 
is the dominant process. At still relatively low laser 
power and low tip voltage, multiphoton photoemission 
over the surface barrier is dominant. As laser power in- 
creases, photoemission from the tip shank contributes. 

This manuscript consists of four main sections. In sec- 
tion II, we explain our experimental setup and our theo- 
retical model. In section HI, we discuss the optical con- 
trol of field-emission sites and the emission mechanism 
based on simulations of local fields on the tip apex and 
laser-induced field-emission microscopy (LFEM) images. 
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In section IV, we discuss the emission processes for vary- 
ing laser power and tip voltage based on experimental 
and calculated results. In the last section, we present 
conclusions and proposals for future experiments. 



II. METHODOLOGY 

A. Experimental setup 

Fig. 1(a) schematically illustrates our experimental 
setup. A tungsten tip is mounted inside a vacuum cham- 
ber (3 • 10~^^ mbar). Laser pulses are generated in a 
Ti:sapphire oscillator (center wavelength: 800 nm; rep- 
etition rate: 76 MHz; pulse width: 55 fs) and intro- 
duced into the vacuum chamber. An aspherical lens (fo- 
cal length: / = 18 mm) is mounted on a holder that is 
movable along the y direction and located just next to 
the tip to focus the laser onto the tip apex; the diam- 
eter of the focused beam is approximately 4 jam (1/e^ 
radius) measured by a method with a razor blade |21| . 
Linearly polarized laser light was used. The polarization 
vector can be changed within the transversal (x, z) plane 
by using a A/2 plate [22]. As shown in the inset, where 
the laser propagates towards the reader's eye as denoted 
by the circled dot, the polarization angle Op is defined 
by the angle between the tip axis and the polarization 
vector. 

The tip can be heated to clean the apex and also nega- 
tively biased for field emission. Although the tip is poly- 
crystalline tungsten, heating to around 2500°C leads to 
the tip apex bein g cr ystallized and oriented towards the 
(Oil) direction (3, l23| . A position sensitive detector with 
a Chevron-type double-channelplate amplifier in front of 
the tip is used to record the emission patterns. The spa- 
tial resolution of field emission microscopy (FEM) is ap- 
proximately 3 nm 0. The tip holder can move alone: 
three linear axes (x, z) and has two rotational axes for 
azimuthal {(p^ around the tip axis) and polar (^, around 
the z axis) angles [IJ. The laser propagates parallel to 
the horizontal y axis within an error of ± 1 degree. is 
set so that the tip axis is orthogonal to the laser propaga- 
tion axis. The orthogonal angle in was determined by 
plotting positions of the tip in (x, y) coordinates while 
keeping the tip apex in the focus of the laser as schemat- 
ically shown in Fig. 1(b). The maximum position in x 
gives the orthogonal angle. Experimental data are shown 
in Fig. 1(c). The plots were taken in 0.5 degree steps. 
The data clearly show the maximum in the x position. 
We defined the corresponding angle as ^ = for conve- 
nience. The precision is estimated to be ± 1 degree. In 
these experiments, the base line of the rectangular detec- 
tor is approximately 20° off from the horizontal {y axis) 
incidence direction, which means that the laser propaga- 
tion axis is inclined by 20° from the horizontal line in 
the observed laser-induced FEM images (see dashed red 
arrow in Fig. 3a). All the measurements were done at 
room temperature. 
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FIG. 1: (color online). Schematic diagram of the experimental 
setup (a). A tungsten tip is mounted inside a vacuum cham- 
ber. Laser pulses are generated outside the vacuum chamber. 
An aspherical lens is located just next to the tip to focus 
the laser onto the tip apex. Emitted electrons are detected 
by a position-sensitive detector in front of the tip. The po- 
larization angle is defined in the inset, where the laser 
beam propagates towards the reader's eye (see text for fur- 
ther description) . (b) shows a schematic diagram defining the 
orthogonal angle between the laser propagation direction and 
the tip axis. The right angle is where x is maximum, (c) 
shows tip positions in (x, y) coordinates for different 0, found 
while the tip apex is kept in the focus of the laser. The angle 
which gives maximum x is defined as ^ = for convenience. 



B. Theoretical model 

Although the field emission is a quantum mechanical 
phenomenon, the interaction between the optical fields 
and the tungsten tip apex can be treated classically by 
solving Maxwell equations. Such an interaction can be 
understood by a mechanistic picture as shown in Fig. 
2(a). When a laser pulse illuminates the metallic tip, 
surface EM waves are excited, which propagate around 
the tip apex. As a result of the interference among the 
excited waves, the optical fields are modulated. To sim- 
ulate the superposition of surface EM waves and the re- 
sulting local field distributions on the tip apex, w e used 
the Multiple Multipole Program (MMP) [25l-[27j. which 
is a highly accurate semi- analytic Maxwell solver, avail- 
able in the package MaX-1 [28| and in the open source 
project OpenMaX [2i,[30j. 

A droplet-like shape was employed as a model tip as 
shown in Fig. 2(b), with a radius of curvature of the 
tip apex of 100 nm, which is a typical value for a clean 
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tungsten tip. Atomic structures were not included in the 
model becasue the tip apex can be regarded as a smooth 
surface on this length scale given by the tip dimensions 
and the wavelength of the laser field. The dielectric func- 
tion e of tungsten at 800 nm was used, i.e. a real part 
Re{e) = 5.2 and an imaginary part Im{e) = 19.4 plj . 
Note that accuracy of the dielectric functions does not 
affect our conclusion, which will be demonstrated in sec- 
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FIG. 2: (color online). Schematic illustration of the exci- 
tation and interference of surface EM waves (a) and of the 
model tip (b). The radius of curvature of the tip apex is 
100 nm, and the length is 700 nm. (c) represents the field 
distribution of the focused beam used in the simulation at a 
certain time. The beam waist is 1 /xm and the wavelength is 
800 nm. Small arrows indicate the field direction, and field 
strength is represented by a linear color scale: highest field 
values are represented in yellow (brightest color). The calcu- 
lated time-averaged field distribution around the model tip is 
shown in (e) for Op — 0° together with a magnified picture in 
the vicinity of the tip apex, (e) shows the time-averaged field 
distribution around a longer model tip together with a mag- 
nified picture in the vicinity of the tip apex; the tip length is 
twice as large as that of (d). (f) represents the time-averaged 
field distributions around the model tip simulated by incident 
laser light with wavelengths of 750 nm, 800 nm and 850 nm. 



tion III B by comparing with local fields on a gold tip. 
A focused laser with a beam waist of 1 /im and a wave- 
length of 800 nm was used as shown in Fig. 2(c). The 
model tip was set so that its apex is at the center of the 
focus. 

By using different droplet sizes it was verified that the 
model tip is long enough so as to mimick the infinite 
length of the real tip; the fields at the truncated side 
of the tip are substantially weaker so that the excited 
surface EM waves propagating arond the whole tip do 
not affect the induced field distribution at the tip apex. 
Figure 2(d) shows the calculated time- averaged field dis- 
tribution around the model tip of Fig. 2(b). Fig. 2(e) 
shows the same calculated field distribution for a longer 
tip with the same radius of curvature of the tip apex of 
100 nm. In both cases, the laser is propagating from left 
to right, where the polarization vector has been chosen 
parallel to the tip axis {Op = 0°). The magnified pictures 
around the tip apex of Fig. 2(d) and 2(e) show that the 
local field distributions of the two are basically the same, 
indicating the length of the model tip in Fig. 2(b) to be 
long enough. 

In the simulations, we only used a center wavelength of 
800 nm, even though the laser pulse has a spectral width 
of A = 25 nm with respect to the center wavelength. Jus- 
tification of the use of only a center wavelength is done 
by simulating the local electric fields with laser light of 
wavelengths of 750 nm and 850 nm. Fig. 2(f) shows 
the time-averaged field distributions around the model 
tip obtained by excitation with these three wavelengths. 
They are almost the same, which indicates that the sub- 
stantial spectral width of the light around 800 nm would 
not affect the position of the maximum local electric fields 
simulated with the wavelength of 800 nm. 



III. OPTICAL CONTROL OF FIELD-EMISSION 
SITES AND EMISSION MECHANISM 

A. Experimental Results 

The field emission pattern from the clean tungsten tip 
apex which orients towards the (Oil) direction is shown 
in Fig. 3(a). The most intense electron emission is ob- 
served around the (310)-type facets, and weaker emission 
from (lll)-type facets. These regions are highlighted by 
green areas with white edges on the schematic front view 
of the tip apex in the inset of Fig. 3(a). The intensity 
map roughly represents a work function map of the tip 
apex: the lower the work function is, the more electrons 
are emitted. The relatively high work functions of (011)- 
and (OOl)-type facets [32[ suppress the field emission from 
those regions. 

The laser-induced FEM (LFEM) image in Fig. 3(b), 
taken with the light polarization oriented parallel to the 
tip axis {Op =0), shows a striking difference in symmetry 
compared to that of the FEM image in Fig. 3(a). Emis- 
sion sites are the same in both cases, but the emission 
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pattern becomes strongly asymmetric with respect to the 
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FIG. 3: (color online). Electron emission patterns for two 
orthogonal azimuthal orientations ((/?) of the tip without laser 
[(f = 0° (a) and cp = 90° (c)], and with laser irradiation 
[(f — 0° (b) and Lp — 90° (d)]. Vup indicates the DC potential 
applied to the tip and Pl indicates the laser power measured 
outside the vacuum chamber. The insets of (a) and (c) show 
the front view of the atomic structure of a tip apex with a cur- 
vature radius of 100 nm, based on a ball model, in which green 
areas with white edges indicate the field emission sites and the 
red arrow indicates the laser propagation direction. The inset 
of (b) shows a schematic side view of the laser-induced field 
emission geometry, in which green vectors indicate intensities 
of electron emission and the white arrow indicates the laser 
propagation direction. A dashed white line denotes a mirror 
symmetry line of the atomic structure in each picture. In 
(c) and (d) specific regions of interest, marked by dashed red 
rectangles, are blown up on the right hand side. 



shadow (right) and exposed (left) sides to the laser inci- 
dence direction. The most intense emission is observed 
on the shadow side as illustrated in the inset of Fig. 3(b). 
Figs. 3(c) and 3(d) give the same comparison for a dif- 
ferent azimuthal orientation of the tip as shown in the 
inset of Fig. 3(c). In the magnified image of Fig. 3(c), 
two emission sites can be identified which are separated 
by approximately 30 nm. The strong emission asymme- 
try is observed even over such short distances, as shown 
in Fig. 3(d). Actually, the laser pulses arrive at an an- 
gle of 20° off the horizontal line in both LFEM images, 
as indicated by the dashed red arrow in the inset of Fig. 
3(a). This oblique incidence slightly affects the symmetry 
with respect to the central horizontal line in the observed 
LFEM images (see below). 

The asymmetry in LFEM images can be controlled fur- 
ther by changing 0. In Fig. 4, the ^-dependence of LFEM 
images at = 0°, 6>p = 0°] is shown, which were taken 
at Vtip = -1500 V and Pl = 20 mW. is varied from 
6> = -12° to 6> = 12° by 4 degree steps. Schematic di- 
agrams for the experimental configuration are shown at 
the top, in which red arrows indicate the laser propaga- 
tion direction. The corresponding FEM images, which 
were taken at Vup = -2200 V, are also shown. As in- 
creases, the asymmetry of the LFEM images becomes 
clearly enhanced. At ^ = 12°, electrons are emitted al- 
most only from right-side emission sites. Among these 9 
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FIG. 4: (color online), ^-dependence of LFEM images at 
[if = 0°, Op = 0°], which were taken at Vup = -1500 V and 
Pl = 20 mW. is varied from 6' = -12° to 6> = 12° by 4° 
steps. Schematic diagrams for the experimental configuration 
are shown at the top, in which red arrows indicate the laser 
propagation direction. The corresponding FEM images are 
also shown below, which were taken at Vup = -2200 V. The 
white dashed lines in the pictures denotes a mirror symmetry 
line of the atomic structure. The total yield S right from right 
side of each image and the total yield from left side Sieft with 
respect to the white dashed line were taken. The ratio of 
S right to Sieft are plotted in the graph. Blue circles are for 
LFEM and black squres are for FEM. 
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FIG. 5: (color online). Comparison of measured and simulated laser-induced FEM (LFEM) images for different light polarization 
angles 6p and for different azimuthal orientations (p of the tip. The leftmost column gives the FEM images without laser 
irradiation for four different azimuthal angles (Vup = -2250 V). For the same azimuthal angles, observed LFEM images are 
shown as a function of polarization angle 9p in 30° steps {Vup ~ -1500 V and Pl = 20 mW). The simulated LFEM images 
from the photo-field emission model, in which Vup and Pl were set as in the corresponding experiments, are shown on the 
right-hand side of the observed LFEM images. The color scale and laser propagation direction are the same as in Fig. 3. 



values, the symmetry of the FEM images changed only 
slightly due to a change of the DC field distribution on 
the tip apex. To distinguish between the contributions 
of DC and laser field distributions to the asymmetry of 
the LFEM images, we evaluated the change in symme- 
try quantitatively. The total yield Sright from the right 
side of each image and the total yield from the left side 
Sieft were obtained from each image with respect to the 
white dashed line. Then the ratios of Sright to Sieft are 
plotted in the graph: high values indicate large asym- 
metries. The asymmetry is clearly enhanced in LFEM 
with respect to FEM, which indicates that the laser fields 
mainly contribute to enhance the asymmetry for higher 
angles 0. 

We also found experimentally a strong dependence of 
the electron emission patterns on the laser polarization 
direction and on the azimuthal tip orientation. Fig. 4 
shows the LFEM patterns for different values of in 
30° steps, and for four different azimuthal orientations 
(p of the tip. The corresponding FEM images are also 
shown in the left-most column; they show simply the 
azimuthal rotation of the low work function facets around 
the tip axis. Throughout the whole image series, the 
emission sites do not change, but intensities are strongly 
modulated resulting in highly asymmetric features. For 
instance, for [cp = 0° ^ 6p = 0°] the intense emission sites 
are located on the right-hand (shadow) side of the tip, 
but for [(p = 0° ^ 6p = 60°] the emission sites on the left- 
hand side become dominant. LFEM images recorded for 
Op = 180 (not shown) are exactly the same as those for 
Op = 0°, and all the LFEM images are well reproducible. 



B. Simulations of Local fields 

When a laser pulse illuminates the metallic tip, surface 
EM waves are excited, which propagate around the tip 



apex as illustrated schematically in Fig. 2(a). Due to 
the resulting interference pattern, the electric fields show 
an asymmetric distribution over the tip apex, depending 
also on the laser polarization. Fig. 6(a) shows the time 
evolution of laser fields at 800 nm wavelength over a cross 
section of the model tip while propagating through the 
tip apex from left to right, where the polarization vector 
has been chosen parallel to the tip axis (6>p = 0°). It can 
be seen that a surface EM wave is propagating around the 
tip apex indicated by white arrows and enhanced at the 
tip apex. The calculated time-averaged field distribution 
around the tip apex is shown in Fig. 6(b). The field 
distribution is clearly asymmetric with respect to the tip 
axis, with a maximum on the shadow side of the tip. The 
field enhancement factor of the maximum point is 2.5 
with respect to the maximum field value of the incident 
laser, and 1.7 for the counterpart of the maximum point 
on the side exposed to the laser. This is consistent with 
our observations in Fig. 3(b) and 3(d) where the field 
emission is enhanced on the shadow side. 

Additionally, we would like to note that a similar asym- 
metric distribution can also be seen even for a metal with 
a dielectric function which is largely different from that 
of tungsten. For example, we performed a simulation for 
gold using a real part Re{e) = — 24 and an imaginary 
part Im{e) = 1.5 [33|. Fig. 6(c) shows the time- averaged 
field distribution on the gold tip apex. The field distribu- 
tion shows a similar asymmetry as for tungsten. It should 
also be mentioned that surface EM waves are classified in 
terms of the dielectric functions of the interacting mate- 
rial [33 though some authours do not distinguish. If the 
real part of the dielectric function is negative, then the 
surface EM waves are proper surface plasmon polaritons. 
On the other hand, if the real part is positive and the 
imaginary part is large, the term Zenneck waves is more 
appropriate. The dielectric functions of tungsten and 
gold between 700 nm and 900 nm are plotted by black 
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FIG. 6: (color online). Time evolution of laser fields over 
a cross section of the model tip while propagating through 
the tip apex from left to right (a). The polarization vector 
is parallel to the tip axis {6p = 0°). Small black arrows 
indicate the field direction, and field strength is represented 
by a linear color scale: highest field values are represented in 
yellow (brightest color). The time-averaged field distribution 
for tungsten and gold tips are shown in (b) and (c) where the 
model tip of Fig. 2(b) is employed for both. The dielectric 
function of tungsten and gold for the wavelengths between 
700 nm and 900 nm are plotted by black dots in (d) and the 
values at 800 nm are highlighted by red circles, (e) shows the 
time-averaged field distributions around the tungsten tip for 
three different polar angles: 6* = -12°, 6' = 0°, and 6' = 12°. 
In (f) the time- averaged field distributions are given in a front 
view of the model tip for different polarization directions 6p 
{0 — 0°). The laser propagation direction is indicated by red 
arrows, and is the same as in the experiment. 



dots in Fig. 6(d), where the values at 800 nm are high- 
lighted by red circles. From Fig. 6(d), strictly speaking, 
the excited surface EM waves on tungsten are Zenneck 
waves and those on gold are surface plasmon polar it ions. 
Figs. 6(b) and (c) also indicate that different kinds of 
surface EM waves do not show substantial difference in 
the resulting field distribution. 

The asymmetric local field distribution can be con- 
trolled by changing the polar angle and the laser po- 
larization angle Op. Fig. 6(e) shows time- averaged field 
distributions on the tungsten tip apex with different laser 
incidence directions relative to the polar orientation of 
the tip apex. As increases, the asymmetry becomes 
stronger. This is consistent with our observation in Fig. 



4 where the most asymmetric emission is observed at 
9 = 12°. Fig. 6(f) shows, in a front view, time-averaged 
field distribution maps from the white dashed line region 
of the model tip in Fig. 6(b). This area corresponds 
roughly to the observed area in our experiments. The 
red arrows indicate the laser propagation direction which 
has been set to the same as in our experimental situation. 
The field distribution changes strongly depending on the 
polarization angle. While the maximum field is located 
directly on the shadow side of the tip for Op = 0°, the 
maximum moves towards the lower side of the graphs in 
concert with the polarization vector for Op = 30° and 
6p = 60°, and reappears on the upper side for Op = 120° 
and Op = 150°. For Op = 90° the polarization vector 
is perpendicular to the tip axis and produces two sym- 
metric field lobes away from the tip apex. In general the 
observed LFEM images show the same intensity modu- 
lations (Fig. 5): each LFEM pattern at 6>p = 30° and 
60° shows pronounced emission at the lower side of the 
image, while each LFEM pattern at 6>p = 120° and 150° 
has maximum emission at the upper side of the image. 



Simulations of LFEM by the photo-field 
emission model 



From the calculated local fields, we further simulated 
the LFEM images by considering the photo-field emission 
mechanism. The current density jcaic of field emission 
can be described in the Fowler-Nordheim the ory based 
on the free-electron model as follows 0, [H, [sH, [36| , 



Jcalc 



em 
2^ 



Wa 



D{W,^,F)f{E)dWdE. (1) 



Here, e is the electron charge and m the electron mass, 
—Wa is the effective constant potential energy inside the 
metal, W is the normal energy with respect to the sur- 
face, and E is the total energy. Important factors are 
D{W, <!>, F) and f{E). D{W, <l>, F) is the probability that 
an electron with the normal energy W penetrates the sur- 
face barrier. It depends exponentially on the triangular- 
shaped potential barrier above W, as indicated by the 
cross-hatched area in Fig. 7(a) where field emission with 
energy W occurs. The cross-hatched area is determined 
by the work function <l> and the electric field F just out- 
side the surface. f{E) is an electron distribution func- 
tion. In the case of field emission we have F = F^c^ 
where Fp,c is the applied DC electric field, and f{E) is 
the Fermi-Dirac distribution at 300 K as shown in Fig. 
7(a). In the photo-field emission model, F still equals 
Fdc, but the electron distribution is strongly modified 
by the electron-hole pair excitations due to single- and 
multi-photon absorption, resulting in a nonequilibrium 
distribution characterized by a steplike profile, as illus- 
trated in Fig. 7(c) [isl,!!^- For example, one-photon ab- 
sorption creates a step of height Si from Ep to Ep^hv by 
exciting electrons from occupied states between Ep — hv 
and Ep. Absorption of a second photon creates a step 
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FIG. 7: (color online). A schematic diagram of field emission 
from a Fermi-Dirac distribution (a), where an electron with a 
normal energy W is emitted. The surface barrier above W is 
shown by a cross-hatched area, (b) shows the logic diagram 
for obtaining work function and absolute DC field maps, (c) 
and (d) show schematic diagrams of photo-field emission from 
a nonequilibrium electron distribution and optical field emis- 
sion from a Fermi-Dirac distribution, respectively. 



of height ^2 from Ej? + hu to Ef + 2hiy^ where ^2 ^ Sf. 
We included absorption of up to four photons. The step 
height Si is proportional to the light intensity /. In the 
vicinity of the tip we have / ex F^^^^^ where Fiaser is 
the enhanced optical electric field that varies over the tip 
apex as illustrated in Fig. 6(f) (ssj . 

There are three adjustable parameters in our calcu- 
lations of jcaic' Fdc^ and ^i. They are all func- 
tions of position on the tip apex. ^ and Fdc maps on 
the tip apex were obtained from the measured FEM im- 
ages by following the logic diagram shown in Fig. 7(b). 
The measured FEM images, which were symmetrized to 
have the ideal two-fold symmetry, represent the current 
density jexp as a function of position on the tip apex. 



TABLE I: Table of work functions of tungsten for several 
faces. An error in experimental values can be as much as 
zb 0.3 eV 0. 



Facet type 


(Oil) 


(001) 


(031) 


Simulation 


4.9 eV 


4.6 eV 


4.45 eV 


Experiment 


5.25 eV" 


4.63 eV^ 


4.35 eV^ 
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FIG. 8: (color online). The calculated DC field distribution 
around the model tip is shown in (a). A wire- like shape was 
employed for the simulation of relative DC fields. The color 
scale is the same as in Fig. 6. The relative DC field distribu- 
tion at the tip apex of (a) is shown as a function of angle Oc, 
which is defined in (a). The obtained work function profile 
along a (001)- (Oil)- (010) curve is shown in (c) as a function 
of Oc. 



because the electrons follow closely the field lines from 
the tip apex to the position sensitive detector. In prac- 
tice we assumed a radius of curvature of the tip apex of 
100 nm and used the FEM image at (/:) = 45° shown in 
Fig. 5. Second, a relative DC field Fdc. relative distribu- 
tion was generated by MaX-1. We used a more wire-like 
tip shape for this purpose as shown in Fig. 8(a), and a 
grounded plate was set 1 cm away from the tip, which is 
close to that in our experimental setup. The simulated 
Fdc. relative is showu in Fig. 8(b), which is normalized 
by the value at tip apex. Going away from the tip apex 
Foc.reiative dccreascs. A scaling factor a is introduced, 
which determines Fdc by Fdc = ol • Fdc. relative- We 
then obtained the <l> map by inserting Fdc into Eq. (1) 
and postulating jexp — jcaic = 0. The resulting <l> map 
was compared to known values for several surface facets 
of tungsten. The scaling factor a was changed and the 
procedure was iterated until reasonable work functions 
were obtained. Thus, a full <l> map and absolute values 
for Fdc were determined. 

A line profile of the resulting <l> map along the (001)- 
(Oll)-(OIO) curve is shown in Fig. 8(c). The work func- 
tion has local curve maxima at the (Oil)- and (001)- type 
facets and local minima at (310)- type facets, which is in 
line with the observed field emission intensity distribu- 
tion seen in Fig. 3(a); the higher the work function, the 
lower the intensity. The resulting ^ values are summa- 
rized for several facets and compared with known experi- 
mental values in Table 1. They are in fair agreement with 
each other. A field strength Fdc of 2.15 V/nm results 
at the tip apex center for the FEM image taken with 
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Vtip = -2250 V, which is a typical value for FEM. The 
LFEM experiments were carried out with a reduced tip 
voltage Vtip ^ —1500 V. Therefore we used a down-scaled 
value of 1.43 V/nm in the LFEM simulations. Note that 
the uncertainty in the <l> values is not important for our 
conclusions which will be discussed below: we have also 
checked the whole discussion in this section with a differ- 
ent work function map using 4.6 eV, 4.32 eV, 4.20 eV for 
(Oil), (001) and (310) facets, respectively, but the main 
outcome does not change. 

Substituting the obtained <l> and F distribution maps 
into Eq. (1), and using a nonequilibrium electron distri- 
bution f{E)^ the absolute values of over the tip apex 
were determined by fitting the measured total current 
from the (310) facet on the right-hand side of the LFEM 
image in Fig. 3(b). The resulting maximum value for 
was 1.6 • 10~^. By substituting all the adjusted pa- 
rameters into Eq. (1), we could simulate all the LFEM 
images. The calculated current densities on the tip apex 
were projected to the flat screen by following the static 
field lines. The simulated images can now directly be 
compared to the experimental images (Fig. 5): they are 
in excellent agreement in every detail. This comparison 
clearly demonstrates that the observed strongly asym- 
metric features originate from the modulation of the local 
photo-fields. 



D. Simulations of LFEM by the optical field 
emission model 



We also simulated the LFEM images for the optical 
field emission process and compared the resulting in- 
tensity distributions to those of the photo-field emission 
model. In the optical field emission model schematically 
shown in Fig. 7(d), the Fermi-Dirac distribution is not 
modified, but instead the electric field F in Eq. (1) is 
expressed as F = F^c + ^laser where F^^^^^ is the nor- 
mal component of Fiaser at each point of the tip surface. 
The absolute values for Ff^^^^ on the tip apex were de- 
termined in the same way as described above for Si . The 
resulting maximum value for F^^^^ was 0.71 V/nm. 

Fig. 9(a) shows the LFEM image for the optical field 
emission process together with experimentally obtained 
LFEM and the simulated LFEM image based on the 
photo-field emission model at [(^ = 0°, 6>p = 0°]. The op- 
tical field emission model results in an even more strongly 
asymmetric pattern as compared to the photo-field emis- 
sion model. This contrasts with the experimental data. 
Fig. 9(b) shows line profiles extracted from the observed 
FEM and LFEM images, and from the corresponding 
simulations for both models, which are all normalized 
to their maximum value. The measured LFEM profile 
clearly shows the asymmetric feature observed in the re- 
gions B and D. The photo-field emission model catches 
this asymmetry much more quantitatively than the op- 
tical field emission model, as can be best seen in region 
D. Therefore, the local fields in our experiment are still 



(a) LFEM 



Photo-field emission Optical field emission 




The number of pixels 

FIG. 9: (color online), (a) experimentally obtained LFEM 
image and simulated LFEM image at = 0°, = 0°] for 
both photo-field emission and optical field emission models, 
(b) shows line profiles extracted from the observed FEM im- 
ages (red line with squares), LFEM images (blue line with 
circles) , and from LFEM images simulated by the photo- field 
emission model (green solid line) and the optical field emission 
model (black dashed line) at = 0°, = 0°]. The whole 
scanned line corresponds to the unfolded rectangle indicated 
by the dashed blue line in the FEM figures above, and the 
corresponding sides are indicated by white arrows. Each line 
profile has been normalized by the maximum value. 



weak enough such that the photo-field emission process 
is the dominant one. 



IV. VOLTAGE- AND POWER- DEPENDENCE 
OF EMISSION PROCESSES 

A. Lower laser power 

In this section, we futher discuss the details of elec- 
tron emission processes and their dependence on laser 
power and tip voltage by investigating electron emission 
patterns and Fowler- Nordheim plots, and by simulating 
electron energy distributions. Figure 10 shows the de- 
pendence of LFEM images at [(/:? = 0°, 6>p = 0°] on the 
average laser power Pl and the tip voltage Vup applied 
to the tip. Throughout the whole seriese of pictures, the 
intensity of each image was normalized by the maximum 
intensity. Normally, total yields decrease as either laser 
power or tip voltage decrease. As can be seen, the left- 
right asymmetry is present in all images below a laser 
power of 60 mW except when = mW (FEM). Nev- 
ertheless, the images show a trend: the outlines of emis- 
sion facets become diffuse in the lower tip voltage region, 
see e.g. the 20 mW row surrounded by a white dashed 
line. In these experiments, electron emission is consid- 
ered to be a concerted action of photoemission and field 
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FIG. 10: (color online). Laser power and tip voltage dependence of the electron emission patterns at = 0°, = 0°]. On the 
vertical axis laser power, on the horizontal axis tip voltage are plotted. The inset shows the electron emission patterns where 
the laser beam is displaced from the tip apex downwards by a distance of 1 /xm and 16 /im, taken at the 90 mW laser power. 
The definition of distance d is also shown on the right-hand side of the inset. The time- averaged field distribution around the 
tip when the laser beam is displaced by a distance of 0.5 /xm is also shown in the inset, where the longer model tip shown in 
Fig. 2(e) was used. The green dashed circles highlight the left-side electron-emission sites (see section IV B). 



emission. In the case of field emission, the emitted elec- 
trons strongly feel the work function corrugation on the 
nanometer scale, which generates a sharp contrasts at the 
border of each emission facet. In the case of photoemis- 
sion, the excited electrons encounter a much narrower 
surface barrier and appear thus to be less sensitive to 
the work function corrugation, hence showing a smeared 
contour of the emission sites. At lower tip voltage, mul- 
tiphoton processes will be enhanced since field emission 
is suppressed. Eventually, photoemission from 3PPE or 
4PPE will contribute, with energies above the vacuum 



level. Therefore, the outline of each emission facet be- 
comes diffuse in the lower tip voltage region. This is 
confirmed by simulations of energy distribution curves 
and emission patterns. 

We have simulated the energy distribution of field- 
emitted electrons with the parameters representing the 
emission site where maximum intensity can be seen in 
the simulated image in Fig. 5 at [(/p = 0°, 6>p = 0°]: the 
work function at this point is 4.45 eV, the DC field is 1.32 
V/nm and is 1.6-10"^, which should correspond to the 
conditions where the most intense point can be seen in 
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FIG. 11: (color online). Simulated electron energy distribu- 
tion curves (a). The spectrum at the top of (a) shows the sim- 
ulated electron energy distribution of field-emitted electrons 
with the parameters representing the point where maximum 
laser intensity can be observed for Lp — 0°, Qp — 0°. The 
parameters are the following: the work function is 4.45 eV, a 
DC field of 1.32 V/nm corresponds to a tip voltage of -1500 V 
and S\ is 1.6- 10~^ for 20 mW laser power. The energy distri- 
bution for lower DC fields but the same laser power are also 
shown in (a). The vacuum level Evac is defined as eV, and 
the Fermi level Ep is 4.45 eV below E^ac The energies cor- 
responding to one-, two- and three-photon excitations from 
Ef (IPPE, 2PPE and 3PPE) are also indicated by vertical 
dashed lines. The simulated LFEM images at corresponding 
tip voltage and laser power are shown in (b). 



the image at 20 mW laser power and -1500 tip voltage in 
Fig. 10. Figure 11(a) shows the corresponding simulated 
energy spectrum at the top, and underneath, spectra for 
various lower tip voltages. We find that field emission 
from two-photon processes is strongly dominant at the 
higher DC fields. For lower DC voltage, the calculated 
energy distributions clearly show that field-emission from 
the 2PPE line is suppressed and photoemission from the 
3PPE line is enhanced. The simulated LFEM images at 
corresponding tip voltage are shown in Fig. 11(b). As the 
DC voltages decrease, the outlines of emission sites be- 
come diffuse due to the fact that photoemission processes 
become dominant. This is in line with the experimental 
observations described above. 

Experimental Fowler-Nordheim (FN) plots give further 
suppports for the suggested emission mechanism. Fig. 
12 shows FN plots of FEM and LFEM, where the elec- 
tron count rate i devided by l^^^ is displayed on a log 
scale versus the inverse of Vup. The count rate was taken 
by integrating the electron emission from the right-hand 
(310) type facet in Fig. 10. According to the FN the- 
ory, the linearity of such plots indicates that the electrons 
are emitted through field emission processes Q . The FN 
plots of FEM data in Fig. 12 clearly show linear be- 
havior, and linearity can be seen also for LFEM at 10 
mW and 20 mW, which are shown together with approx- 
imated exponential functions (black solid lines). The FN 
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FIG. 12: (color online). Fowler-Nordheim (FN) plots for FEM 
and LFEM {Pl = 10 mW, 20 mW, 30 mW, 40 mW, 50 mW). 
The vertical axis is signal i over the on a logarithmic 
scale. The signal % is the total yield of electron emission from 
the right hand (310) type facet of each image in Fig. 10. The 
horisontal axis is l/Vtip. 



plots for 20 mW show non-linear behavior at low bias 
voltage. This indicates that photoemission processes be- 
come dominant at low voltage as discussed above. Sim- 
ilar behavior is also observed in the case of higher laser 
power, 30 mW, 40 mW and 50 mW. 

The slope of the straight sections is proportional to 
<|)3/2 From this fact, we can estimate the effec- 

tive barrier height ^LFEM{\^mW) and ^LFEM{2^mW) 
dit Pl =10 mW and 20 mW, respectively, which 
an emitted electron feels in the case of LFEM. First, 
the barrier height ratios were derived from the propor- 
tionality constants: ^LFEM{iOmW)/^FEM = 0.24 and 
^LFEMi20mW)/^FEM = 0.2. Taking the work func- 
tion of (310) type facets of 4.35 eV for ^fem^ thus 

^LFEM{10mW) = 1-05 cV and ^LFEMi20mW) = 0-85 cV 

were obtained. The energy difference between ^lfem 
and ^FEM should corresponds to the energy of emit- 
ted electrons measured from the Fermi level in LFEM. 
Here we obtain ^LFEMiiQmW) - ^fem = 3.3 eV and 



LFEM{20mW) 



FEM 



3.5 eV, which is close to 



the electron energy after two-photon excitation, i.e. 3.1 
eV. These values corroborate the two-photon photo-field 
emission processes, which is consistent with the simu- 
lated electron energy distributions for higher voltage in 
Fig. 11(a). 



B. Higher laser power 

For this discussion, we would like to point out the 
electron emission from the shank side of the tip in the 
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higher laser power range. As in the previous sections, 
the electron emission from the tip apex is dominant in 
the emission pattern at low laser power because of the 
local field enhancement at the tip apex even though the 
surface area exposed to the laser beam is much larger for 
the shank than for the apex. At the higher laser power, 
however, the electron emission from the shank side be- 
comes noticiable because of the nonlinear dependence of 
the electron emission intensities on the laser power. 

In the column of Vup = — lOOV of Fig. 10, the left-side 
electron-emission features highlighted by green dashed 
circles becomes suddenly very strong for laser powers ex- 
ceeding 70 mW. At 90 mW, the intensity of left-side elec- 
tron emission is comparable to that on the right-side. It 
remains even when the position of the tip in the beam 
waist is varied. The insets of Fig. 10 show electron emis- 
sion patterns at 90 mW laser power where the laser beam 
is displaced from the tip apex downwards by distance of 
1 fim and 16 /im. In the two images, right-side emission 
sites disappear, but the left-side emission remains, indi- 
cating that it originates from the shank side of the tip. 
Such an electron emission should be dominated by pho- 
toemission over the surface barrier because DC fields on 
the shank side are significantly weak with respect to the 
tip apex. Since the laser pulses arrive at an angle of 20° 
off the horizontal line in both LFEM images, the position 
of the electron emission from the shank is also deviated 
from the horizontal line. In the inset, we also show the 
time-averaged field distribution around the tip when the 
laser beam is displaced downwards from the tip apex by 
a distance of 0.5 /im: the longer model tip shown in Fig. 
2(e) was used. The maximum field can be observed at 
the side exposed to laser, which is consistent with our 
observations. 



V. CONCLUSIONS 

We have observed laser-induced modulations of field 
emission intensity distributions resulting in strong asym- 
metries, which originate from the laser-induced local 
fields on the tip apex. By varying the laser polariza- 



tion and the laser incidence direction relative to both 
azimuthal and polar orientation of the tip apex, we have 
demonstrated the realization of an ultrafast pulsed field- 
emission source with convenient control of nanometer 
sized emission sites. These experimental observations 
are quantitatively reproduced by using simulated local 
fields for the photo-field emission model. We discussed 
the emission processes further and found field-emission 
after two photon photo-excitation to be the dominant 
process in laser-induced field emission. From experimen- 
tal data and simulations, the dependence of the emission 
processes on laser power and tip voltage could be under- 
stood. 

This type of electron source is potentially useful for 
many applications such as time-resolved electron mi- 
croscopy, spatio-temporal spectroscopy, near-field imag- 
ing techniques, surface-enhanced Raman spectroscopy, or 
coherent chemical reaction control. Maybe the most in- 
teresting applications will arise when two laser pulses 
with different polarizations or paths are used for the emis- 
sion of two independent electron beams from two differ- 
ent sites on the tip, spaced only a few tens of nanometers 
apart, and with an adjustable time delay between the two 
electron pulses. Since field emission electron sources pro- 
duce highly coherent electron beams due to their inher- 
ently small source size, comparable to the finite spatial 
extent of electron wave packets inside the source 3, Q , we 
could expect two spatially and temporally coherent elec- 
tron beams to be available within the coherence time. 
This should create new opportunities for addressing fun- 
damental questions in quantum mechanics such as anti- 
correlation of electron waves in vacuum [4l|, or for new 
directions in electron holography (42j . 
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